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The nature of the superconducting transition in highly underdoped thick films of La 2 -a;Sra;Cu 04 
{x = 0.07 and 0.08) has been investigated using the in-plane transport measurements. The contribu¬ 
tion of superconducting fluctuations to the conductivity in zero magnetic field, or paraconductivity, 
was determined from the magnetoresistance measured in fields applied perpendicular to the Cu 02 
planes. Both the temperature dependence of the paraconductivity above the transition and the 
nonlinear current-voltage (/ — V) characteristics measured across it, exhibit the main signatures of 
the Berezinskii-Kosterlitz-Thouless (BKT) transition. The quantitative comparison of the super¬ 
fluid stiffness, extracted from the 1 — V data, with the renormalization-group results for the BKT 
theory, reveals a large value of the vortex-core energy. This finding is confirmed by the analysis 
of the paraconductivity obtained using different methods. The results strongly suggest that the 
characteristic energy scale controlling the BKT behavior in this layered system corresponds to the 
superfluid stiffness of a few layers. 

PACS numbers: 74.40.-n, 74.72.Gh, 74.62.En 


I. INTRODUCTION 

One of the most intriguing phenomena in condensed 
matter systems is the occurrence of the so-called 
Berezinskii-Kosterlitz-Thoulessi^— (BKT) transition in 
two-dimensional (2D) superfluid systems. The main in¬ 
gredients of the BKT physics were described originally 
within the context of the two-dimensional XY model, 
which is an effective model for the collective phase of the 
superfluid order parameter.i^ Here logarithmically inter¬ 
acting vortex-like topological excitations drive the tran¬ 
sition from the superfluid state, where they are bound 
together in vortex-antivortex (V-AV) pairs, to the metal¬ 
lic one, where single vortex excitations proliferate. This 
mechanism leads in principle to several peculiar signa¬ 
tures in the physical observables, such as the univer¬ 
sal and discontinuous jump^ of the superfluid density at 
Tbkt, the observation of which in ^He films^ was the first 
experimental proof of the existence of a BKT transition. 
Afterwards, interest in BKT physics was triggered mainly 
by the possibility to observe it in superconducting (SC) 
systems that can be considered to be in the 2D limit. 
On very general grounds, this occurs for systems with 
low superfluid stiffness Jg, defined as the energy scale 
associated to the areal density of superfluid electrons: 

_ h^nsdsKT _ dsKT /..X 

“ 16^e2 A2 ’ 

where rig , m denote the superfluid density and mass of 
the carriers, respectively, A is the magnetic penetration 
depth and dsKT denotes a transverse length scale over 
which the system can be seen as effectively 2D. The pos¬ 
sibility to see BKT physics is connected to a low value of 
dsifr/A^: indeed, despite the presence of screening su¬ 
percurrents, the interaction between vortices remains log¬ 


arithmic when the Pearl screening length A = 2}?/dsKT 
overcomes the system sizei^ In addition, since the dis¬ 
tance between Tbkt and the ordinary BCS tempera¬ 
ture Tc scales as [Tc — Tbkt)/Tc oc Tbkt/ Js, a clear 
BKT regime can only be identified when Jg gets reduced. 
In films of conventional superconductors, these condi¬ 
tions are usually realized when the film thickness d is 
reduced. In those cases, by identifying dBKT with d, typ¬ 
ical BKT signatures have been observed^— by means of 
different experimental probes. The universal jump of the 
superfluid density has been seen either via direct mea¬ 
surements of the inverse penetration depth^riiii^ or via 
a discontinuous jump of the exponent of the nonlinear 
I — V characteristicsiS At the same time, the vortex pro¬ 
liferation above Tbkt has been identified^ d"^d^ from an 
exponential divergence of the correlation length above 
Tbkt, which leads to a peculiar paraconductivity above 
the transition Jdi 

An alternative route for the observation of BKT 
physics is presented by bulk layered systems, in which 
the magnetic-field distribution of a vortex differs drasti¬ 
cally from the monopole-like Pearl solution in uniform 
films the presence of other superconducting layers 
squeezes the field of a pancake vortex into a narrow strip 
of size A along the c axis. This in turn implies that the 
logarithmic dependence of the interaction potential be¬ 
tween two vortices placed in the same layer persists up 
to all length scales, as in a neutral superfluid, making 
in principle the stack of uncoupled layers the best pos¬ 
sible system to observe a true BKT transition, with the 
2D unit in Eq. ([T]) corresponding to each isolated plane. 
In the presence of Josephson coupling between layers, 
the upper cut-off for the logarithmic interaction between 
vortices becomes^d^ Aj ~ ^o/\/J±/J\\, where is the 
zero-temperature in-plane coherence length, and J|| _l are 
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the in-plane and out-of-plane superfluid stiffness, respec¬ 
tively. If the interlayer coupling is weak, i.e. Jl/J|| 1, 
this length scale is large enough to allow for a BKT-like 
description of the vortex-antivortex interaction, indepen¬ 
dent of the film thickness d. In practice, even if the finite- 
size effect due to A,/ leads to a rounding of the discontin¬ 
uous jump in Js, the analysis of anisotropic 3D XY-\ike 
modeli^Ti^i shows that the unbinding of vortex-antivortex 
pairs in each plane is still the mechanism driving the tran¬ 
sition, in analogy with the purely 2D case. Therefore, in 
a weakly coupled, layered superconductor, one expects 
to observe a BKT-like transition at a 3D transition tem¬ 
perature that is slightly higher than the BKT transition 
temperature of a single layer of an equivalent uncoupled 
system. 

Such a description is expected to be appropriate for 
underdoped samples of cuprate superconductors, which 
are highly anisotropic, layered materials. Here one also 
finds that the superfluid stiffness is suppressed by the 
proximity to the Mott insulatormaking the separa¬ 
tion between Tbkt and Tc large, while avoiding the addi¬ 
tional consequences of an increase of the disorder level, as 
it occurs in films of conventional superconductors when 
the thickness is reduced. According to this argument, 
in bulk samples of underdoped cuprates one should be 
able to identify BKT signatures assuming that the fun¬ 
damental 2D unit is represented by isolated Cu02 layers, 
i.e. the transverse length scale dsKT in Eq. 0 would 
coincide with the interlayer distance dc^ as pointed out 
in the seminal work by Emery and Kivelson^ However, 
it has been recently shown^^ that this picture is some¬ 
how too simplified, since one should also account for the 
nontrivial role of the vortex-core energy p, which is the 
energetic cost needed to create the vortex at the small¬ 
est length scale ^o- Indeed, even if the layers are weakly 
coupled, what matters for the vortex proliferation is the 
competition at large distances between the effective vor¬ 
tex fugacity and the effective Josephson coupling. As 
a consequence, when p is large, the Josephson coupling 
between layers can prevent the vortex unbinding, mov¬ 
ing the BKT transition away from the value expected for 
each isolated layer, resulting in an effective dimension 
dsKT larger than dc. 

So far, the experimental situation in cuprate super¬ 
conductors has been controversial. For example, the di¬ 
rect measurements of the inverse penetration depth have 
shown that, in the YBa 2 Cu 307 _a; family, no BKT jump 
is observed even in strongly-underdoped thick films21i2^ 
or crystalsi2S A BKT-like superfluid-density jump is 
only seen in few-unit-cell thick films of YBa 2 Cu 307 _a, 
(Ref. [sO) or Bi9Sr9CaCuT0«^.,.fRef. [^. but even in this 
case, as the samples get underdoped, the effective dsKT 
seems to cross over to the sample thickness and the 
superfluid-density jump gets smeared out. While this 
can be explained indeed by an increase of the vortex- 
core energy with underdopingj^^ one should notice that 
the simultaneous appearance of an anomalously large dis¬ 
sipative response suggests that spurious finite-frequency 


effects can also be present, as emphasized recently in the 
analysis of thin films of NbN These spurious effects are 
instead absent in the dc measurements of the I — V ex¬ 
ponent that suggested a BKT-like jump very near Tc in 
cuprate samplesHowever, this measurement allows 
one to extract directly the effective 2D areal stiffness o, 
i.e. the combination dsKT/^'^-, so dsKT can be deter¬ 
mined only if is known by measurements in similar sam¬ 
ples. Finally, the analysis of the paraconductivity, i.e. of 
the SC fluctuations above Tc, also raises some questions 
on the occurrence or not of a BKT transition. Indeed, 
on one hand, the SC fluctuations have been proved to 
have a str ong 2D character in several cuprate families 
{e.g. Refs. I39l - l43h with the typical 2D unit being identi¬ 
fied as the distance between the Cu02 layers dc- On the 
other hand, these are ordinary Gaussian (amplitude and 
phase) fluctuations, with a BKT regime that, if present, 
is restricted to a small range of temperatures near Tc in 
the most underdoped samples 1 ^ 2 *^ 

In the present work, we address the issue of the identi¬ 
fication of the scale d bkt in cuprate superconductors by 
making a simultaneous analysis of the BKT signatures 
both below and above Tbkt in two highly underdoped 
samples of La 2 -a;Sr 3 ;Cu 04 . We first extract the paracon¬ 
ductivity above Tbkt (Sec. IllBl) . and then determine 
the temperature dependence of the anomalous 2D expo¬ 
nent of the I — V characteristics across it (Sec. |HC|. 
In Sec. IIII Al the direct comparison of the experimental 
I — V data with the renormalization-group results for 
the BKT theory allows us to extract a large value of the 
vortex-core energy p, consistent with that obtained from 
the analysis of paraconductivity in Sec. Ill Bl According 
to earlier theoretical work)22i^ the large value of /i ob¬ 
tained in our study corresponds to dBKT — (2 — 3)dc- 
Furthermore, this value of the vortex-core energy can 
be used to reduce considerably the fitting parameters in 
the well-known Halperin-Nelson formula^ for the para¬ 
conductivity above Tbkt, spanning both the BKT and 
Aslamazov-Larkin^— (AL) regimes of the SC fluctua¬ 
tions. This analysis (Sec. IHIBI) confirms that the effec¬ 
tive length scale dBKT is a few times larger than dc, in 
agreement with the expectatio n^^i^^ for a layered weakly- 
coupled system with a large vortex-core energy. Our 
study clarifies how different transverse length scales en¬ 
ter in the analysis of the SC fluctuations above and below 
Tbkt, solving the apparent contradiction between pre¬ 
vious measurements. 


II. EXPERIMENTS 
A. Samples and measurement techniques 

The samples were La 2 -a;Sra,Cu 04 (LSCO) films with 
the nominal doping x = 0.07 and x = 0.08. They 
were patterned into standard Hall bars with the length 
L = 2.0 mm and the width W = 0.3 mm; the distance 
between voltage contacts was I.Ol mm. The films were 
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75 unit cells (150 Cu 02 layers) thick (d 1000 A) and 
grown by molecular beam epitaxy. The films and samples 
were described in detail elsewhere.— The samples become 
superconductors below the temperature Th^q^x), defined 
as the temperature at which the in-plane resistance R be¬ 
comes zero. The measured Tr=o were (3.9 ± 0.1) K and 
(9.7 ± 0.3) K for samples x = 0.07 and 0.08, respectively. 

The in-plane sample resistance and magnetoresistance 
were measured in ^He cryostats (base T Ri 0.25 K) with 
a standard four-probe ac method (~13-16 Hz) in the 
Ohmic regime, using either SR7265 lock-in amplifiers or 
a LR-700 resistance bridge. The magnetic fields H up to 
18 T were applied perpendicular to Cu 02 planes {H \\ c 
axis) and swept at constant temperatures. The sweep 
rates of 0.02-1 T/min were low enough to avoid the heat¬ 
ing of the sample from eddy currents. 

The current-voltage (I-V) measurements were carried 
out at constant temperatures (T) in 77 = 0 using ^He and 
variable-temperature insert (base T r; 1.3 K) cryostats. 
DC square pulses provided by a Keithley 6221 current 
source were applied to the samples, while a Keithley 
2182A nanovoltmeter measured the voltage response. 
Each data point on the I-V curve was found by averaging 
measurements with positive and negative pulse polarities. 
Such a four-point dc method^ avoids possible effects of 
parasitic capacitances {e.g. from the sample contacts) 
and obviates Joule heating, while retaining the increased 
sensitivity of a finite-frequency technique and eliminat¬ 
ing the effects of thermal electromotive forces. Current 
excitations between 50 nA and 1 mA were typically used, 
depending on the film doping and temperature. 

The addition of current noise to a device with an in¬ 
trinsic nonlinear behavior can create an Ohmic response 
at low currents^ and, in particular, it can create Ohmic 
behavior even below Tbkt- Therefore, for the I-V mea¬ 
surements, filtering was provided at room temperature 
by a 1.75 nF low-pass tt filter in series with a 1 kO re¬ 
sistor on each lead to the sample. The tt filters and the 
resistors were encased in a shielded box attached to the 
top of the cryostat probe. This filter box provided a 5 dB 
(60 dB) noise reduction at 10 MHz (1 GHz), which en¬ 
abled the observation of nonlinear I-V behavior at low 
excitations amid masking current noise. 


B. High-field magnetoresistance measurements and 
superconducting fluctuations 

By approaching the superconducting transition from 
above, it is in principle possible to identify the BKT tran¬ 
sition from the temperature dependence of the contribu¬ 
tion of superconducting fluctuations (SCFs) to conduc¬ 
tivity (or “paraconductivity”), AascF{T) = p{T)~^ — 
Pn{T)~^, where p{T) and Pn(T) are the measured and 
normal-state resistivity, respectively. In cuprates, the de¬ 
termination of Pn(T) has been somewhat ambiguous and 
controversial (see, e.g., Ref. and references therein). 
We emphasize, however, that the precise determination 


of (finite) is not crucial for the extraction of AascF 
in the regime of interest, very near the BKT transition 
where the contribution of SCFs diverges (Eq. ([3]) below). 
On the other hand, it may introduce considerable er¬ 
rors into the values of AascF far from it— This issue is 
demonstrated and discussed further in Sec. mm 

In this Section, we adopt a method that uses trans¬ 
verse (H II c) magnetoresistance measurements to deter¬ 
mine the extent of SCFs. In particular, above a suffi¬ 
ciently high magnetic field i7'(r), SCFs are completely 
suppressed (i.e. they become unobservable in the ex¬ 
periment) and the normal state is fully restored. In the 
normal state, the magnetoresistance of cuprates increases 
as at low field a^^i^^i^^~ — {oJcT <C 1 , where oJc is the 
cyclotron frequency and r is the scattering time), similar 
to the classical orbital effect in conventional metals— 


Pnjll) - pnjO) 

Pn{0) 


oc 77^, 


( 2 ) 


Therefore, the values of 77' can be found from the down¬ 
ward deviations from such quadratic dependence that 
arise from SCFs when 77 < 77 ' 43 ^^^^ 
tribution to the conductivity can be determined then as 
AascFiT, 77) = p{T, H)-^ - p^{T, H)-\ where p{T, 77) 
is the measured resistivity and pn{T,H) is obtained by 
extrapolating the region of 77^ magnetoresistance ob¬ 
served at high enough 77 and T. The advantages of this 
method^i^ over some of the earlier ones {e.g. Refs. 
and are that it does not rely on any assumptions 
about the T dependence of pn, and it makes it possible to 
determine both the paraconductivity AascF{T, H = 0 ) 
and the SCF contribution to conductivity in the presence 
of magnetic field. 

Figure [T]( a) shows representative p(77) curves (77 || c) 
obtained on the x = 0.07 LSCO sample. The condition 
for the weak-held limit is satished in the entire regime of 
interest as uJcT r 0.5 at 18 T and 5 K, where it reaches 
its maximum value. By tracking the gradual evolution 
of the magnetoresistance curves measured at different T 
[Fig- [11(b)], from the high-T region where the 77^ depen¬ 
dence is unambiguous, to lower T where SCFs are more 
pronounced, we were able to determine the values of the 
onset fields 77'(T) (see Appendix lAl for a more detailed 
discussion). 

Figure UKa) inset shows 77'(r), determined from 
Fig. [T](b) for the x = 0.07 sample and fitted by a sim¬ 
ple quadratic expression II{.{T) = 77'(0)[1 — (r/r 2 )^], 
similar to earlier studies— In zero field, SCFs 
become observable below T 2 = 29 K. In Sec. IlHBl we 
show explicitly that the exact determination of and 
thus the determinations of 77'(T) and T 2 , does not affect 
our conclusions. Hereafter we focus only on the zero-field 
behavior. 

Figure UKa) shows that AascFiH = 0)/cr„(77 = 0), 
where (t„ = 1/pn, increases by several orders of mag¬ 
nitude as temperature is reduced towards Tfi=o R! 4 K, 
reminiscent of the exponential divergence expected at the 
BKT transition. Indeed, in 2D the paraconductivity can 
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FIG. 1. (Color online) x = 0.07 LSCO film, (a) Resistivity 
vs. transverse magnetic field {H || c) up to 18 T for differ¬ 
ent T, as shown. The highest T data are also shown in the 
inset for clarity, (b) Magnetoresistance data from (a) plot¬ 
ted vs H^. Symbols (black diamonds) show H'^iT), the fields 
above which SCFs are fully suppressed and the depen¬ 
dence of the normal-state resistivity pn is observed. Dashed 
lines are linear hts representing the contributions from normal 
state transport, i.e. they correspond to [p{H) — p(0)]/p(0) = 
[pn(0) - p(0)]/p(0) -I- [pniO) / p{0)]atranaH^ ■ 
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FIG. 2. (Color online) x = 0.07 LSCO film, (a) Symbols 
show /S.gscf{T, H = 0)/(t„( 77 = 0) vs T, as determined from 
Fig. [T] The solid line is a ht to Eq. (|4]) with Tbkt = 3.8 K 
and fitting parameters A = 17.9, h = 2.9; the dashed line 
corresponds to Tbkt = 4.0 K, A = 29, b = 3.3. Inset: The 
onset field H'i.{T) below which SCFs become observable. The 
dashed line is a fit 77' = 77'(0)[1 — (r/T 2 )^], with po77'(0) = 
15 T and T 2 = 29 K. (b) p vs T in 77 = 0. Arrows point at 
Tbkt= Tr^o and Tc; Tc was estimated as shown in the inset. 
Inset: The temperature at the inflection point of the p(T) 
curve, where dp/dT has a maximum, is taken as an estimate 
of Tc in the calculation of a from Eq. ©• 


always be expressed as 


^<ySCFl<^n 


'aT)Y 

. ^0. 


(3) 


where ^(T) is the SC correlation length, whose tempera¬ 
ture dependence depends on the nature of the SC fluctu¬ 
ations. The usual Aslamazov-Larkin^— (AL) paracon- 
ductivity describes the fluctuating Cooper pairs above 
the mean-field temperature Tc, and leads to a power-law 
divergence of the coherence length ~ (T — Tc)~^. In 
contrast, within BKT theory, ~ l/np measures 

the inverse density np oi free vortices above Tbkt, and 
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diverges exponentially as T —>■ Tbkt- An interpolation 
formula between these two regimes was first proposed by 
Halperin and Nelson^ 

where t = {T — Tbkt)/Tbkt, and A and b are numerical 
constants. More recently, a renormalization-group (RG) 
study^ of the BKT transition showed that parameter b 
is strictly connected to two relevant physical quantities: 

b ~ 2a\/t/, a = ijl/ijlxy, (5) 


where tc is the distance between the mean-field and BKT 
critical temperatures 


_ Tc — Tbkt 
Tbkt 


( 6 ) 


while a is the vortex-core energy /r expressed in units of 
the conventional value /rxv that it assumes in the XY 
model (see also Eq. (IT6l) below). According to Eq. (|3l), 
the exponential BKT behavior is limited to the range of 
temperatures t <C tc, while above it, one recovers the 
usual AL paraconductivity. 

The paraconductivity shown in Fig. [2Ka) has been fit¬ 
ted to Eq. Q by taking Tbkt = Tb=o = (3.9 ± 0.1) K 
[Fig [2Kb)]. Surprisingly, it is possible to get a good fit 
to the data even up to very high temperatures ^ 20 K 
with reasonable values of A and b [e.g. dashed line in 
Fig.|2Ka)]. However, within the error for Tbkt, the lower- 
T data up to ^ 10 K are described better with the fitting 
parameters in the range A = 13—20 and b = 2.5—3.0 [e.g. 
solid line with A = 17.9 and b = 2.9 in Fig.|2Ka)]. Assum¬ 
ing that Tc ^ 7 K, i.e. of the order of the temperature 
where dp/dT has a maximum [Fig |2Kb) inset], Eq. ([S]) 
then yields enhanced values of the vortex-core energy, 
p/pLxY — 1.4— 1.7, consistent with previous workj ^^i^^ 

The above analysis of the SCFs above a SC transition, 
which occurs at Tbkt = Tb.=o, suggests the presence of 
a BKT fluctuation regime at Tbkt < T < Tc ~ 7 K, 
followed by a crossover to the AL regime at T > Tc. It is 
worth noting that the crossover to the AL regime gives 
some indication on the transverse length scale controlling 
the Gaussian fluctuations in the sample. Indeed, when 
tc, Eq. ([31) reduces to 


sc F 


4 &^ Tbkt 
T — Tbkt 


— KbKT 


Tc 

T-Tc’ 


(7) 


where, on the r.h.s., kbkt = APTbkt /[A^Tc) and we 
replaced T — Tbkt with cs T — Tc, which is correct when 
T is sufficiently larger than Tc so that the difference be¬ 
tween Tc and Tbkt can be neglected. We note that, in 
films of conventional superconductors usually tc is at 
most of order 0.1, so the crossover from the pure BKT be¬ 
havior to the AL one occurs for relatively small reduced 
temperatures t. In our samples, tc is as large as 0.7, 
so the asymptotic AL behavior (I7|) is reached at higher 


temperatures. On the other hand, since the SCFs regime 
extends up to reduced temperatures as large as t ~ 3 — 4 
[Fig.|2Ka)[, there is still a large temperature regime where 
the approximation ([7|) is valid. This expression has to be 
compared with the usual AL formula^^^— that gives 

AuAB _ PnjdxL Tc _ Tc 

“ IQRc T-Tc^'^^^T-Tc’ 

where Rc = h/e^ = 4.1 kO. By mapping the expres¬ 
sions ([7|) and (|8|), we can see that the high-T limit of the 
interpolating HN formula also hxes the prefactor kal 
that controls the strength of the AL fluctuations in the 
Gaussian regime at T > Tc. The latter one depends 
in turn on the transverse length scale dxL that identi¬ 
fies the 2D unit for AL fluctuations [see Eq. ([5|)[. By 
using the estimates of b, A given above, we obtain that 
Kbkt = Aj^Tbkt/{A^T c) ~ 0.1. Thus, from the mea¬ 
sured ~ 1 mil cm and by matching kbkt and kal , we 
conclude that dAL is of the same order as the interlayer 
distance dc, in full agreement with previous work in the 
literature4I^ In other words, as far as the Cooper-pair 
fluctuations are concerned, the fluctuation regime dis¬ 
plays marked 2D character with decoupled layers, consis¬ 
tent with the standard expectation for a weakly-coupled 
layered superconductor 42 On the other hand, the BKT 
paraconductivity does not allow us to extract any pre¬ 
cise information on the scale dBKT controlling the vor¬ 
tex physics below Tbkt- To address this issue, and to 
confirm the fit based on the paraconductivity data ex¬ 
tracted from the high-held magnetoresistance measure¬ 
ments [Fig. HJa)], we analyze the I — V characteristics, 
whose behavior is, in fact, one of the key signatures of 
the BKT transition. 


C. Current-voltage characteristics and superfluid 
stiffness 

The most famous hallmark of the BKT transition is 
observed by approaching Tbkt from below. In particu¬ 
lar, the superfluid stiffness Jg, dehned in Eq. exhibits 
the so-called universal jump at the transition, i.e. 

Js{Tbkt) ~ ~Tbkt, Js{T^j^j,) = 0. (9) 

TT 

Here the T dependence of Js(T) includes both the quasi¬ 
particle excitations, which would drive Jg continuously to 
zero at Tc, and vortex-like phase fluctuations, which are 
instead responsible for the discontinuous jump ([9]). The 
latter directly influences the behavior of the exponent a 
in the I — V characteristics: 

a(T) = + 1. (10) 

The superlinear behavior in Eq. m is due to the ability 
of a sufficiently large current to unbind vortex-antivortex 
pairs. From Eq. ([21), it follows then that a should jump 
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from a = 3 at T = to a = 1 at T = 

Below Tbkti the exponent a is expected to increase with 
decreasing T since the superfluid density increases. 

The voltage-current characteristics are shown in 
Figs. Efa) and [3)Jb) on a log-log scale for the x = 0.07 
and X = 0.08 films, respectively. The power-law behav¬ 
ior V oc is observed at all T in the low-current 

limit. In that regime, the V{I) dependence is thought 
to arise from the thermally dissociated vortex-antivortex 
pairs for T > Tbxt ^-nd from current-induced dissocia¬ 
tion for T < Tbkt- At the highest currents in Figs.[3l^a) 
and[3Jb), heating effects become important. The temper¬ 
ature dependent exponents a(T) were determined as the 
slopes of the linear fits of the data at the lowest currents 
[Figs.[3fa) and[3fb)]. We note that, due to a large value 
of tc, the fitting range, both in current and in tempera¬ 
ture, is much wider than usual, i.e. compared to systems 
that are clearly 2D, such as interfaces^ and films.— The 
values of a{T) are presented in Fig.[3ljc) for both samples. 
A steep change of a from its Ohmic value (a = 1) at high 
T to large values > 3 is indeed observed with decreasing 
T. In particular, a(T) in the x = 0.08 sample exhibits a 
jump-like behavior as expected, but the a(T) dependence 
is smoother in a more highly underdoped sample. Never¬ 
theless, a reaches 3 at T = (3.6±0.1) K and (9.4±0.1) K 
for samples x = 0.07 and x = 0.08, respectively, close to 
their Tfj=o values and consistent with the assumption in 
Sec. Ill B I that = Tbkt- 

Even though the Tbkt values will be determined more 
precisely in Sec. IIII Al by the theoretical analysis that 
takes into account the smearing of the BKT jump by the 
presence of inhomogeneities, we can estimate the order 
of magnitude of Js{Tbkt) from the temperature where 
a(T) = 3 using Eq. (jH]). In the x = 0.07 sample, for 
example, we have Js{Tbkt ~ 3.6 K) Ri 2.3 K. Using this 
value and Eq. m expressed as^ 


J^K] = 0.62 


dBKT [A] 
A2[/im2] ’ 


( 11 ) 


we find that, if the effective transverse length scale 
dBKT coincides with the film thickness (« 10^ A), 
^{Tbkt) ~ 16 /iin, while for dBKT — dc 7 A we obtain 
^{Tbkt) ~ 1.4 pm. Based on the doping and tempera¬ 
ture dependences of the penetration depth measured in 
similar LSCO filmS)^ we estimate that X{Tbkt) does not 
exceed a value of 2-3 pm for our x = 0.07 sample. There¬ 
fore, we find much better agreement between the results 
of our I—V measurements and penetration depth studies 
by assuming that the effective sample thickness is some¬ 
what larger than the interlayer spacing, but not as large 
as the whole thickness of the sample. As we shall see be¬ 
low, this conclusion is confirmed by a detailed comparison 
between Js(T) extracted from the a(T) exponent and the 
theoretical prediction for the BKT behavior, when the 
non-trivial role of the vortex-core energy is taken into 
account. 

Einally, we remark that, in our samples, we do not 
expect to observe the Ohmic response in the I — V char- 



/(A) 
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FIG. 3. (a), (b) Voltage-current characteristics on a log-log 
scale for x = 0.07 and x = 0.08 samples, respectively, at differ¬ 
ent T, as shown. At the lowest excitations, V oc I°'V) ^ where 
the solid lines are linear fits with the slopes corresponding to 
a{T). In both panels, the dashed lines with slopes a = 1 and 
a = 3 gnide the eye. (c) a[T) for both samples. The dashed 
line a = 3 crosses the data at (3.6 ± 0.1) K and (9.4 ± 0.1) K 
for samples x = 0.07 and x = 0.08, respectively. 
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acteristics caused by finite-size effectsi^^— Indeed, it is 
known that the Ac I — V curves probe the contribution 
of dissociated vortex-antivortex pairs separated by a dis¬ 
tance r* = 2'kJscW/^qI. Therefore, at small currents, 
which probe r* larger than the sample width {W < L), 
the free vortices will dominate the resistance and the 
I — V characteristics will be Ohmic. On the other hand, 
the nonlinear behavior m of the I — V characteristics 
can only be seen when r* < W, i.e. above a threshold 

currentiiSa /*, 

_ 2Js7rc ^ AksTsKT c 

“ $0 “ $0 ■ 


the presence or not of free-vortex excitations, described 
by the large-distance behavior of the vortex fugacity. The 
physical superfluid stiffness Jg is then obtained by the 
numerical solution of the RG equations at large distances 
(see Appendix [Bl for more details). 

Apart from the starting value of which can 

be determined by comparison with the data far from 
Tbkt, the second relevant energy scale in the problem is 
the ratio fj./ Js- Here we take it as a free parameter, to 
be determined by the fit to the experimental data. This 
has to be contrasted to the usual AT-model description 
of the BKT transition, where fi/Jg is constrained to the 
value 


By using the above estimate Tbkt ~ 3.6 K for the x = 
0.07 sample, one gets I* ~ 2.68 x 10“®(A/K) TsifriK] ~ 
1 X 10”^ A. In the presence of inhomogeneous domains of 
size L' < L, the threshold current I* is expected^ to in¬ 
crease with respect to the estimate m- However, since 
the homogeneous value m we found is considerably 
smaller than the currents at which the measurements are 
performed, finite-size effects are not expected to man¬ 
ifest themselves in our experiment. Indeed, Figs. Ha) 
and Ht>) show that, below Tbkt, the crossover from 
the nonlinear behavior (flUl) back to the Ohmic one is not 
observed even at the lowest measured current. 


III. THEORETICAL ANALYSIS OF THE DATA 


A. Superfluid stiffness 


We extract from Eq. m the temperature depen¬ 
dence of the superfluid stiffness Js{T), which we analyze 
along the lines of the approach discussed earlier for both 
conventiona l and cuprate superconductors . ^^,32 
In Eq. (fTUl) . the temperature dependence of the super¬ 
fluid stiffness Js(T) is due to both quasiparticle excita¬ 
tions, which induce a BCS-like suppression of 
at all temperatures up to Tc, and vortex-like excitations, 
which become relevant near Tbkt < 2c- Since our I — V 
measurements are rather close to Tc, we can assume for 
jBCSfrp-^ a linear behavior, 

J^t^S(T) = Jo ■ (13) 


The effect of vortices is taken into account by solving 
the BKT renormalization-group equations, whose rele¬ 
vant variables are 


K = 


7rJ^<^^(r) 

T 


g = 27re 


(14) 

(15) 


where g is called the vortex fugacity (/3 = l/fegT). Here 
determines the value of K at the shortest length 
scale of the problem, i.e. the SC coherence length ^q, 
while the large-distance behavior will be determined by 


TT 

fJ-XY = (16) 

In general, the value of g/ Jg determines the tempera¬ 
ture scale where significant deviations of Jg from the 
BCS temperature dependence J^^^{T) start to become 
visible. Indeed, even though free vortices only start to 
proliferate at Tbkt, if a significant density of vortex- 
antivortex pairs already exists below Tbkt, it can renor¬ 
malize {i.e. suppress) the large-distance superfluid stiff¬ 
ness Jg{T) with respect to its BCS behavior counterpart 
much before the BKT transition. In thin films of con¬ 
ventional superconductors it has been shown that this 
is the caseJ^i^ Here fi/Jg ~ 1, as expected in ordinary 
BCS superconductors, and the measured Jg{T) deviates 
from the BCS behavior significantly before the univer¬ 
sal jump (H occurs. In contrast, it has been argued^^i^ 
that, in cuprate superconductors, g/Jg can even exceed 
the (large) value ~ 4.9 in Eq. (II6II . where Jg is now the 
stiffness of a single layer {i.e. with Abkt = dc in Eq. ([1])). 
As we shall see below, this has relevant consequences for 
the determination of the effective transverse scale Abkt 
for the BKT transition in a bulk material or in a thick 
film, as it is in our case. 

A second effect to be taken into account in the analysis 
of the experiments is the presence of inhomogeneity of the 
local SC properties, which have been clearly shown to be 
relevant in underdoped cuprates by means of STM analy¬ 
sis of underdoped samplesi ^^d^ Here we modeli^i^i^ the 
presence of inhomogeneities by assuming that the local 
BKT critical temperature has a finite distribution about 
the most probable value, represented by the curve labeled 
“Homogeneous” in Fig. 01 The main effect of the inho¬ 
mogeneity is then to smear out the universal jump ([5]). 
the effect being larger for a wider probability distribu¬ 
tion of the local Tbkt values. More details are given in 
Appendix IbI 

The results for the two samples x = 0.07 and x = 0.08 
are shown in Fig. 01 and the fitting parameters are sum¬ 
marized in Table HI The BCS temperature dependence 
m, shown in the figure with a dotted line, reproduces 
the data below the BKT transition very well, in partic¬ 
ular for the X = 0.08 sample where more experimental 
points are available. Here the sample inhomogeneity is 
very small (width of the distribution S/Jq = 0.01; see also 










FIG. 4. Temperature dependence of the superfluid stiffness 
for the X = 0.08 (a) and x = 0.07 (b) samples: Comparison 
between the experiment and the theory, as described in the 
main text. The fitting parameters are listed in Table H] 


TABLE I. Fitting parameters for Fig. [I] 


doping 

Jo (K) 

T. (K) 

Tbkt (K) 

/^/ Js 

5/Jo 

tc 

^theo 

0.07 

6.5 

6.5 

4 

6.3 

0.1 

0.625 

2.02 

0.08 

41 

11.3 

9.7 

7 

0.01 

0.16 

1.15 


Appendix [B]) and, accordingly, the homogeneous and in¬ 
homogeneous curves almost coincide, with a sharp down¬ 
turn of Js(T) near Tbkt- Tbkt is defined here as the 
transition temperature for the homogeneous curve, which 
represents the most probable transition temperature for 
the sample. We note that, since we also included the ef¬ 
fects of the finite size of the system, which lead to some 
rounding of the Js{T) jump before Tbkt-, even in the 
homogeneous case we do not observe a strictly discon¬ 
tinuous jump as in Eq. (jOj, but Jg vanishes continuously 
over a temperature range of a few mK. For the x = 0.07 
sample, the inhomogeneity is larger {5/ Jq = 0.1), as ex¬ 
pected for a more underdoped sample, and this leads in 


particular to a longer superfluid tail above Tbkt- In 
both samples, we extract a large value of the vortex-core 
energy, i.e. ^/Js = 6 — 7 or ^xy ~ 1-4. As explained 
above, this implies that the deviations of Js{T) from the 
BCS curve only occur near Tbkt- As a consequence, 
Tbkt can be very well estimated by using the universal 
relation ([HI) with Js{Tbkt) replaced by J^^^iTgKT)^ 
i.e. 

== ^ 3^, (17) 

TT TrJo 

which is in very good agreement with the tc values listed 
in Table U obtained by the RG results. It is apparent 
that the large separation between Tc and Tbkt in our 
samples is due to the presence of two concomitant ef¬ 
fects in underdoped cuprate films: (i) the large mean-field 
critical temperature and (ii) the low superfluid stiffness, 
proportional to Jq in Eq. 03, due to correlations 
This has to be contrasted to conventional superconduc¬ 
tors, where the BKT regime can only become visible 
when Jq is suppressed by strong disorder, which also 
brings along unavoidable spurious effects connected to 
the inhomogeneityjHi^iSi In addition, in systems like 
NbN, it has been shown that fx/Js — 1, so the devia¬ 
tions of Js{T) from the BCS behavior occur much before 
the intersection with the BKT line,— making the ap¬ 
proximate estimate (fTTl) much less reliable. 

Our finding of the large value of /r is an important re¬ 
sult, since it confirms previous theoretical analysi o^^i^^ 
in cuprates, and it allows us to understand the estimated 
value of d-BKT ^ dc in our film, as discussed in Sec. Ill Cl 
Since the measurements of a(T) only access the areal 
superfluid stiffness (HD and thus do not allow for a sep¬ 
arate determination of cIbkt and A, the comparison of 
the experimental data and the theory shown in Fig. |4] 
has been done for the BKT transition in the pure 2D 
case. On the other hand, we also know that our films 
are comprised of ^ 10^ layers, with a weak interlayer 
Josephson coupling between them. In this case, it has 
been proven by previous theoretical wor k^^i^^ that, when 
fJ-/fJ-XY > Ij the BKT transition Tbkt moves away from 
the value expected for a single, isolated layer In 

particular, according to the analysis of Refs. 1231 and 1^. 
for the value of /r found above, one could expect that 
Tbkt about 30% larger than corresponding to 

dBKT — (2 — 3) dc. Indeed, by assuming dBKT — ‘2-dc for 
the X = 0.07 sample, for example, one can easily estimate 
Tbkt = ~ 1.3using the r.h.s. of Eq. (IT71) . 

The value dBKT — (2 — 3) dc is consistent with the esti¬ 
mate based on the comparison to the penetration-depth 
measurements discussed in Sec. Ecl 

B. Paraconductivity 

We note also that the value of fj. extracted from the 
behavior of Js{T) is consistent with that obtained in 
Sec. Ill BI from the analysis of the paraconductivity above 
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TbkTi even though the fits presented there do not in¬ 
clude the effect of SC inhomogeneities. Indeed, we can 
show that for our samples the inhomogeneity has a rela¬ 
tively minor effect on the determination of the parame¬ 
ters entering the paraconductivity fit. To show this, we 
analyze the paraconductivity above Tbkt by refining the 
analysis of Sec. Ill Bi with the inclusion of inhomogeneity. 

We can describe the measured resistivity as 


A 

Rn 


A ^crsCF 
V O-n 


-1 



( 18 ) 

In order to compare the theoretical predictions to a 
larger number of data points, in AascpiT) = p(T)~^ — 
Pn{T)~^, we approximate pn = l/on with a constant, 
zero-field value measured at T ^ Tbkt- Even though 
this procedure is less accurate far from Tbkt, this is 
not relevant for the discussion of the effects near Tbkt, 
where the SCFs contribution diverges. This is exem¬ 
plified in Fig. [SJ where paraconductivity obtained using 
this method is compared to that extracted from mea¬ 
surements in high magnetic fields (Sec. Ill Bl) . Indeed, the 
exact determination of pn becomes crucial only far from 
Tc, where Agscf j becomes comparable to the differ¬ 
ence 10%) between its values obtained using those 
two methods. 

For the x = 0.07 sample, we choose = i?(r = 
20 K), where the SCFs contribution to conductivity is 
only a few per cent [Fig. [5]) a)]. To account for the inho¬ 
mogeneity, we follow the procedure discussed in Ref. H . 
and outlined in Appendix |BJ We use the same distri¬ 
bution of local critical temperatures extracted from the 
analysis of Js{T) to generate a distribution of local re¬ 
sistivity values Ri/Rn described by Eq. (1181) with the 
same local values T^kt^^^c computed above. Thus, only 
A, b in Eq. (1181) are the fitting parameters. The global 
resistance of the sample is then determined by the cor¬ 
responding random-resistor-network problem by means 
of the effective-medium approximation. Once again, to 
elucidate the role of inhomogeneity, we compare the re¬ 
sults for the homogeneous and inhomogeneous case. The 
“Homogenous” curve in Fig. [S] refers to the paraconduc¬ 
tivity of a system with a single Tbkt and T^ realization, 
corresponding to the most probable value in the sample. 
Thus, this is the paraconductivity expected for a homo¬ 
geneous system whose superfluid stiffness below Tbkt 
is described by the “Homogeneous” Js{T) curve in Fig. 
Hb). 

The results are shown in Fig.[S]for the parameters A = 
14 and b = 2.55, which are in good agreement with the 
results of the analysis in Sec. imi The differences are 
due to the effect of the inhomogeneity, which is known^i 
to affect the slope of R{T) above the transition. More 
importantly, 6 = 2.55 is very close to the theoretical value 
btheo = 2 calculated from Eqs. (ED, i.e. 


btheo — 9 T 

Js 


(19) 



FIG. 5. Comparison between the R{T)/Rn experimental 
data and the theoretical prediction obtained in the homoge¬ 
nous or inhomogeneous case. In the latter case, R{T) is ob¬ 
tained as solution of a random-resistor-network problem in the 
effective-medium approximation, as explained in Appendix m 
Solid diamonds: Paraconductivity was extracted from mea¬ 
surements in high magnetic fields fSec. Ill Bl) . 


by using the values of p/Js and tc extracted from the 
analysis of the I — V characteristics below Tbkt, and 
listed in Table H] The fit accurately reproduces the data 
up to T ~ 10 K, which is an extremely large range of 
SCFs, similar to Fig. [IJa). However, the BKT fluctua¬ 
tion regime only extends up to Ri 6.5 K and, above it, 
ordinary AL-like Gaussian fluctuations are at play. Fi¬ 
nally, we note that some deviations start to occur above 
T ~ 10 K. As we discuss in Appendix |b 1 this effect has 
already been observed in several families of cuprates^i^ 
and it can be interpreted as a signature of a pseudogap 
state above 


IV. DISCUSSION 

The analysis carried out in the previous Sections 
clearly demonstrates the occurrence of a BKT-like tran¬ 
sition in our LSCO samples. This is confirmed both by 
the analysis of the paraconductivity above Tbkt and by 
the analysis of the superfluid stiffness below Tbkt, as 
extracted from the I — V measurements. Even though 
our films are thick, in the sense that d is much larger 
than the SC coherence length, the possibility to see BKT 
physics is guaranteed by the layered nature of the sys¬ 
tem. As we discussed in Sec. HI a weakly-coupled layered 
superconductor is an ideal candidate for observing the 
BKT physics, since a layered structure ensures the best 
screening of the charged supercurrents. Indeed, in this 
case the interaction between vortices in each plane is log¬ 
arithmic up to a scale Aj ~ \fJ±TJ\\ that grows as 
the stiffness anisotropy increases. While the behavior of 
J_l/J|| as a function of doping in the LSCO family has 
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not been systematically explored, in other cuprates it has 
been shown to decrease significantly with underdoping)^^ 
along with a general suppression^ZiSl of J|| due to cor¬ 
relation effects i^^'^^ Under these conditions, one could 
expect to identify signatures reminiscent of the typical 
2D BKT physics, such as an almost discontinuous sup¬ 
pression of the superfluid stiffness, even in a layered bulk 
samplef^i2£i2^i2i On general grounds, the starting point 
of this reasoning is that, as demonstrated within several 
modelsi^Ti^i the physics of a layered superconductor with 
a very weak interlayer coupling closely approaches that 
of an isolated 2D system. Indeed, even though the tran¬ 
sition will ultimately have a 3D character, the 3D critical 
region is extremely reduced for weak interlayer Josephson 
coupling,— >21 and it could even be masked in the experi¬ 
ments due to finite-size effects or inhomogeneities of the 
type discussed in this manuscript. 

Since in the BKT picture there exists a universal re¬ 
lation (ini) between the transition temperature and the 
smallest superfluid stiffness beyond which vortex unbind¬ 
ing occurs, the idea that each layer is isolated can lead 
to the naive expectation that Tbkt is controlled by the 
stiffness of each isolated layer, i.e. the value © with 
dsKT = dc- However, as predicted theoretically^ this 
simple picture should be in part revised when the role of 
the vortex-core energy /r, controlling the vortex fugacity 
g ~ is taken into account. Indeed, in a layered 

BKT model the transition temperature is not controlled 
by the “bare” {i.e. short-distance) values of J±/J\\ and of 
the vortex density g, but by their large-distance behavior. 
Both energy scales grow at large distances, with opposite 
consequences: the increasing of J±/J\\ tries to keep the 
system superconducting, while the increase of g implies 
that vortices would like to proliferate making the system 
non-superconducting. While at some temperature g will 
ultimately win, the counter-action of the interlayer cou¬ 
pling can move Tbkt away from the temperature scale 
connected to the single-layer stiffness. Thus, the effective 
stiffness to be used in Eq. ([S]) has to be computed from 
the definition o with a transverse length scale dBKT 
somewhat larger than dc. In particular, as g, increases, 
the transition temperature moves farther away from the 
single-layer temperature scale— The large value of the 
vortex-core energy obtained in our measurements sug¬ 
gests that, in strongly-underdoped LSCO samples, the 
relevant length scale controlling the BKT transition in¬ 
volves a few coupled layers, i.e. dBKT ^ dc- This 
conclusion is in agreement with the estimate based on 
the measured superfluid stiffness o, i.e. a combination 
dBKT /^^1 and the comparison to A measured in similar 
films. 

These findings, based on the analysis of the super- 
fluid stiffness below Tbkt, are confirmed by the anal¬ 
ysis of the paraconductivity above it. In particular, we 
have shown that the SC fluctuations above Tbkt exhibit 
a BKT character near the transition, and then evolve 
into the ordinary Aslamazov-Larkin-type behavior ex¬ 
pected for Gaussian (amplitude and phase) fluctuations. 


We fitted the data with the well-known Halperin-Nelson 
formula)2i which interpolates between the two regimes, 
by constraining the fitting parameters according to the 
theoretical expectations for them— This procedure not 
only provides a consistency check of the validity of the 
BKT analysis, but it also allows us to confirm the esti¬ 
mate of the vortex-core energy extracted by the analysis 
of superfluid density. In agreement with previous findings 
in bulk cuprates most of the fluctuation regime 

is dominated by Gaussian fluctuations with a marked 
2D character, where the characteristic 2D unit is rep¬ 
resented by a single layer, i.e. dAL = dc. It is worth 
stressing that this result is not in contradiction with the 
finding dBKT ^ dc for the BKT behavior. Indeed, in 
the case of Gaussian fluctuations, the dimensionality of 
the fluctuations is controlled only by the band-parameter 
anisotropy, i.e. the ratio t_L/t|| between interlayer and in¬ 
tralayer hopping, respectively. When this ratio is small, 
as it is in cuprates, one can see 2D fluctuations over a 
a broad temperature range— The crossover to 3D be¬ 
havior, expected in bulk materials, is here preceded by 
the vortex fluctuations, which drive the system towards a 
2D BKT transition. Even if the transition will ultimately 
have a 3D character, we do not identify the crossover to 
3D fluctuations. This is consistent with the fact that 
the 3D critical regime, especially above the transition)^! 
is extremely reduced in a weakly-coupled layered system 
and, in addition, it gets masked by inhomogeneous effects 
that are mostly relevant at the transition. 


V . CONCLUSIONS 

We have presented measurements of the in-plane trans¬ 
port properties of two strongly underdoped thick films of 
La 2 -xSr 2 ;Cu 04 . Our results have (i) established the oc¬ 
currence of a BKT-like transition and (ii) identified the 
typical transverse length scale that defines the equivalent 
two-dimensional unit controlling the BKT signatures in 
this layered system. 

The most striking signature of a vortex-driven phase 
transition emerges from the superfluid stiffness Jg, ex¬ 
tracted from the exponent of the nonlinear I — V char¬ 
acteristics across Tbkt- In both samples, we observe a 
rapid downturn of Jg reminiscent of the well-known uni¬ 
versal jump expected in a 2D superconductor. A quanti¬ 
tative comparison with the theoretical predictions, which 
also include the effect of some unavoidable degree of in¬ 
homogeneity in the samples, strongly suggests a large 
energetic cost to create the vortex cores in the SG state. 
As a consequence, even though the interlayer Josephson 
coupling is weak, the vortex-pair unbinding occurs at a 
temperature larger than the one where each isolated layer 
would undergo the BKT transition— In other words, the 
characteristic energy scale controlling the BKT proper¬ 
ties corresponds to the superfluid stiffness of a few lay¬ 
ers. These results are confirmed by the analysis of the 
paraconductivity above Tbkt- Thanks to the few-K dis- 
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tance between the BKT (Tbkt) and mean-field (Tc) crit¬ 
ical temperatures, we can clearly see that an initial BKT 
regime of fluctuations crosses over to an extended regime 
of 2D Aslamazov-Larkin-type Gaussian fluctuations. 

As we remarked above, the advantage of using highly 
underdoped thick films is that the intrinsically low value 
of the superfluid stiffness, due to the proximity to the 
Mott-insulating phase^^i^^ allows us to achieve a large 
separation between Tbkt and Tc without simultane¬ 
ously introducing a large disorder-driven inhomogeneity 
of the local SC properties. This has to be contrasted 
with the case of few-unit-cell thick films of cupratesj^i^ 
which are usually much more sensitive to disorder, so 
that the BKT jump of the superfluid stiffness is usu¬ 
ally lost with underdopingi^ We note also that finite- 
frequency probes, such as the two-coil mutual induc¬ 
tance technique used in Ref. 1^ . can be potentially much 
more sensitive to disorder-induced inhomogeneity, as dis¬ 
cussed recently within the context of films of conventional 
superconductors In contrast, the superfluid density ex¬ 
tracted from the I—V characteristics is a purely dc probe, 
and this can explain why we see a relatively sharp BKT 
jump even in our strongly underdoped samples. Whether 
these features are common to other families of cuprates 
is an interesting open question that certainly deserves 
further experimental and theoretical investigation. 


ACKNOWLEDGMENTS 

We thank A. T. Bollinger and I. Bozovic for the sam¬ 
ples. We acknowledge V. Dobrosavljevic for useful dis¬ 
cussions. High-field {H > 9 T) measurements were per¬ 
formed in the National High Magnetic Field Labora¬ 
tory (NHMFL) DC Field Facility. This work was par¬ 
tially supported by NSF Grants No. DMR-0905843, No. 
DMR-1307075, and the NHMFL through the NSF Coop¬ 
erative Agreement No. DMR-1157490 and the State of 
Florida. L.B. acknowledges financial support by MIUR 
under projects FIRB-HybridNanoDev-RBFR1236VV, 
PRIN-RIDEIRON-2012X3YFZ2 and Premiali-2012 AB- 
NANOTECH. 


Appendix A: Magnetoresistance 

The dependence of the magnetoresistance is clearly 
observed at the highest T and H [Fig.[ljb)]. As the tem¬ 
perature is lowered and SCFs become stronger, the re¬ 
gion gets pushed to higher fields and the curvature of the 
p{H) dependence at high H, in the normal state, becomes 
less obvious. The same kind of behavior has been ob¬ 
served in other cuprates, e.g. in YBa2Cu307_a; (Ref. 
and in overdoped^^*^!^ underdoped^ and even non- 
superconducting^i^ LSCO very close to the onset of su¬ 
perconductivity. 

In underdoped LSCO, it is well know n^^i'^'^~ — that the 
resistivity at high H increases with decreasing T {i.e. 


dp/dT < 0), as seen also in Fig. [T](a), reflecting the ten¬ 
dency towards an insulating ground state at high fieldsj^ 
Nevertheless, deviations from Eq. ([2|) still provide a good 
estimate of iL'(r), as discussed below. While the precise 
reason for the applicability of Eq. m remains an open 
problem beyond the scope of this study, we note that, in 
the regime of interest, the system remains in the (poor) 
metallic regime, as kpl ^ 1 {kp - Fermi wave vector, I - 
mean free path), i.e. the resistance per square per Cu02 
layer < hje^. 

The quadratic dependence H'^{T) = i7'(0)[l —(r/r2)^] 
(see Fig. [2Ka) inset) was found also in YBa2Cu307-.r 
(Refs. 1^ and [H) and overdoped LSCO (Ref. giv¬ 
ing us further confidence that the values of are 

reliable. Furthermore, the H = 0 onset temperature for 
SCFs, T2 = 29 K, is consistent with the results from tera¬ 
hertz spectroscopy^ obtained on similar films, and those 
determined from the onset of diamagnetism^ and the 
Nernst effectS^ in LSCO crystals with similar p{T) and 
Tc values. We also find that poHc{0) ~ 15 T is in agree¬ 
ment with the value of the upper critical field obtained 
from specific-heat measurements^^ on LSCO with a sim¬ 
ilar Tc. Therefore, even though the magnetoresistance 
method that we employed to determine iL' and Pn{H) 
has an inherent limitation in accuracy, we conclude that 
both the magnitude and the temperature dependence of 
the onset fields iJ' are fairly consistent with those from 
other types of studies. This consistency check confirms 
further that the observed onset of the magnetoresis¬ 
tance corresponds to the return to the normal state. 


Appendix B: Renormalization-group analysis of the 
BKT transition for an inhomogeneous system 


The BKT RG equations describe the large-distance be¬ 
havior of the dimensionless quantities K and g introduced 
in Eqs. (fHl) - (fT51) above. They are given by^iiidi 


dK 

~dl 

d£ 


-K^9\ 

i2-K)g, 


(Bl) 

(B2) 


where £ = Inr/^o is the rescaled length scale with respect 
to the short-distance cut-off for the problem, represented 
by the SC coherence length ^o- The initial values of K 
and g are determined by the BCS value of the superfluid 
stiffness, Eq. m, which includes only the temperature 
dependence due to quasiparticle excitations. The effect of 
vortices is accounted by the RG flow at large distances, so 
that the physical superfluid stiffness (HD is identified with 
the limiting value of K as one goes to large distances!^ 


Js = 


TK{£ ^ 00) 

TT 


(B3) 


The basic idea of the RG equations is to look at the 
competition at large scales between the superfluid stiff¬ 
ness and the vortex fugacity. When g —>■ 0, it means that 
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single-vortex excitations are ruled out from the system, 
which then remains superconducting. Indeed, as one can 
see from Eqs. (IBlII and (IB2II . when g ^ 0, K goes to a 
constant and then Jg from Eq. (|B3I) is hnite. If instead 
g —>■ c» at large distances, it means that vortices prolif¬ 
erate and drive the transition to the non-SC state, since 
K ^ 0. The large-scale behavior depends on the ini¬ 
tial values of the coupling constants if, g, which in turn 
depend on the temperature. The BKT transition tem¬ 
perature is defined as the highest value of T such that K 
flows to a hnite value, so that Js is hnite. This occurs 
at the hxed point K = 2, g = 0, so that at the transition 
one always has K{i = oo) = 2, corresponding to the uni¬ 
versal relation ([9]) quoted above. By numerically solving 
Eqs. (|Bl(l - (jB2(l at each temperature, while taking Tc and 
g./ Jbcs 8.S free parameters in the initial value, we obtain 
the curve labeled as “Homogeneous” in Fig. |4l with the 
parameters reported in Table HI 

To account for the presence of inhomogeneities, we fol¬ 
low the procedure discussed in previous publications for 
both conventional^^ and cuprate^ superconductors. We 
assume that the BCS superhuid density is described by 
Eq. (|13|) with the initial value Jg randomly distributed 
according to a probability density P{Jq) that we take, 
for instance, as Gaussian: 

P{J^) = ^ exp [-(J* - Jo)V2S^] . (B4) 

In the homogeneous case, the Gaussian distribution has 
zero width and only the value Jg is allowed. In this case, 
one obtains the Js{T) curve labeled as “Homogeneous” 
in Fig. 01 and the corresponding Tc, Tbkt are the ones 
reported in Table HI As we remarked in the text, we also 
add hnite-size effects, by stopping the RG flow at the 
scale imax = A/^g of the system size. As a consequence, 
even for the homogeneous case Js{T) does not display 
a real jump, but an extremely rapid downturn occur¬ 
ring over a few-mK temperature range. In the inhomo¬ 
geneous case, for each Jg value distributed according to 
Eq. (IB4I) . we rescale the corresponding BCS tempera¬ 
tures as Jg/T^ = Jq/Tc a-nd we compute Jl{T) and the 
corresponding BKT temperatures by the numeri¬ 

cal solution of the RG equations (fT^ - dTSll above. After 
obtaining this set of Jl{T) curves, we compute the sam¬ 
ple stiffness as the average one Jav, dehned as 

JaviT)=Y,P(.JomT). (B5) 

i 

When all the stiffness values Jl{T) are different from 
zero, as it is the case at low temperatures, the aver¬ 
age stiffness will be centered around the center of the 
Gaussian distribution (IB4I) . so that it will coincide with 
jBCS^rpy However, by approaching Tbkt dehned by 
the average not all the patches make the tran¬ 

sition at the same temperature, so that the BKT jump is 
rounded and Jav remains hnite above the average Tbkt, 
in agreement with the experiments. In this analysis, we 


then have a second free parameter that is the width d/Jo 
of the Gaussian distribution (IB4I) . However, all four pa¬ 
rameters of the ht (average Jg and Tc, ratio p/Jg and 
5/Js) affect in a rather independent way the shape of 
the overall stiffness. Indeed, Jg and Tc are essentially 
determined by the slope of the stiffness before the BKT 
transition, /i/Jg determines the location of the univer¬ 
sal jump, whose smearing is controlled by 5. Thus, even 
though some hexibility is possible in the values of the pa¬ 
rameters listed in Table HI these variations are expected 
to be within 10% — 20% of the quoted values. 

The inhomogeneity also influences the paraconductiv- 
ity above Trkt- To show this, we proceed in analogy 
with Ref. [ij by mapping the spatial inhomogeneity of 
the sample in a random-resistor-network problem. In 
particular, we associate to each patch of initial stiff¬ 
ness value Jg a normalized resistance pi = Ri/Rn ob¬ 
tained from Eq. m by using the corresponding local 
values of T^ and T^bkt computed as outlined above. 
The overall sample normalized resistance p = R/Rn is 
then calculated in the so-called effective-medium-theory 
(EMT) approximatioufSi where p is the solution of the 
self-consistent equation 

Rip-p,) 

. P + P^ 

Here Pi is the occurrence probability of each resistor, 
which coincides with the distribution function (jB4ll of 
the local Jg value used to compute the corresponding 
Pi{T). The resulting R{T) = pRn is shown in Fig.jSJ and 
it is compared to the one of the homogeneous case, i.e. 
the R{T) curve obtained when only the most-probable 
Jg value of the distribution (IB4I) is realized. 

As we observed in Sec. HU for T > 10 K the ex¬ 
perimental paraconductivity saturates more rapidly than 
what is predicted by the HN interpolating formula. Since 
in this regime we are already exploring Gaussian fluctu¬ 
ations, such a failure is not correlated with the BKT 
character of the fluctuations, but it pertains instead to 
the regime of ordinary Cooper-pairs fluctuations. In¬ 
terestingly, such behavior has been already observed in 
several families of cuprates and it has been inter¬ 
preted theoretically^!^ as an effect of the pseudogap. 
Indeed, by phenomenological modelling of the suppres¬ 
sion in the electronic density of states characteristic of 
a preformed pseudogap, one can reproduce a faster de¬ 
cay of the Cooper-pairs correlation length ^(T) in Eq. 
(01) with respect to the standard AL prediction. Even 
though a detailed analysis of this issue is beyond the 
scope of the present manuscript, we nonetheless observed 
that a similar effect seems to be at play also in the case 
of our sample. To account for it within the HN interpola¬ 
tion scheme, we can for instance multiply the correlation 
length entering the paraconductivity formula by a 
function suppressing it around a temperature T* larger 
than Tc, such as 

^(T) = j sinh A exp(-(T/T*)4). (B7) 
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By introducing this correction factor in each normalized 
resistivity pi appearing in Eq. (IB6I) . we obtain the (ho¬ 
mogeneous and inhomogeneous) curves displayed in Fig. 
ini Here we used T* = 19 K, that is, approximately the 
temperature where magnetoresistance saturates. As one 
can see in Fig.[6l our scheme now gives an excellent agree¬ 
ment with the experimental data up to T ~ 20 K. In this 
high-temperature regime, the deviations of R{T) from 
the magnetic-field extracted paraconductivity (symbols 
in Figs. [S] and [B]) become sizeable, and the expression 
(IB7p reproduces the latter points well. 



FIG. 6. Comparison between the R{T)/Rn experimental data 
and the theoretical prediction obtained in the homogenous or 
inhomogeneous case for the modified HN correlation length 
(iBTll . including phenomenologically the pseudogap effect in 
the AL fluctuations. Solid diamonds: Paraconductivity ex¬ 
tracted from measurements in high magnetic fields ISec. lII Bll . 
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